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PARTICLE TRACKABILITY CONSIDERATIONS FOR LASER 
DOPPLER VELOCIMETRY 
by Allen M. Karchmer 
Lewis Research Center 

SUMMARY 

The technique for calculating single -particle trackability in a fluid flow in terms of 
rms velocities is discussed. The general equation for the motion of a single particle as 
a response to a fluid -flow driving force is presented, and a criterion for greatly simpli- 
fying the governing equation is developed. It is shown that in the case of large particle- 
to fluid-density ratio and low-frequency fluid motion compared to the natural time con- 
stant of the particle response, the governing equation can be reduced to a simple balance 
of inertia force and linear viscous drag (i.e. , Stokes drag). 

INTRODUCTION 

The advent of the laser Doppler velocimeter (LDV) has presented to the experimen- 
ter interested in making velocity measurements in flowing fluids an opportunity to make 
such measurements without the limitations and disadvantages of previous instruments. 
One of the more significant advantages of the LDV is its relative insensitivity to temper- 
ature and density, thus making the separation of compressibility effects from velocity 
signal unnecessary. There are a number of other advantages, but since there is an 
abundance of literature on this subject (ref. 1, e.g.), they are not listed in this report. 

The principle behind velocity measurements in a fluid from Doppler -shifted laser 
light is well established and straightforward. Light from a laser beam is directed into 
a moving fluid. A portion of the beam is scattered from particles moving with the fluid. 
This scattered light is slightly shifted in frequency by the Doppler effect. The magni- 
tude of the frequency shift is then measured and can be directly related to the velocities 
of the scattering centers. The scattering centers (either artifically introduced or na- 
turally present) must satisfy two basic criteria: they must be efficient scatterers of the 
laser radiation and they must be accurately following, or tracking, the fluid flow. Sev- 



eral authors have shown that the scattering criteria require the LDV to be collecting ra- 
diation scattered in the Mie regime (ref. 1, e.g.). Basically, this amounts to having 
the diameter of the scattering center (i.e. , the seeding particle) at least as large as the 
wavelength of the incident light X. For an argon laser, X is 0.488 micrometer (4880 A); 
therefore, the particle diameters must be approximately 1/2 micrometer or larger. 

This requirement puts an approximate lower limit on the size of the seeding parti- 
cles used. So, except for particles whose density is the same as that of the fluid, there 
will always be some minimum tracking error. This error must be estimated for the 
particular flow of interest to obtain an overall estimate of the accuracy of any LDV mea- 
surements. This report presents techniques to enable individuals involved in laser - 
Doppler velocity measurements to estimate the errors involved in their measurements 
because of finite particle response. 


SINGLE-PARTICLE RESPONSE TO FLUID FLOW 
Particle Motion as Response to Driving Force 

If the particle velocity is considered the output, or response, to an input driving 
force of a turbulent fluid velocity, the transfer function is defined through the relation 

v p (jcu) = G(j w)v f (j u>) (1) 


where Vp(ja>) is the Fourier transform of the particle and v f (ju>) is the Fourier trans- 
form of the fluid velocity. All symbols are defined in the appendix. The energy spec- 
trum of the particle velocity is | v (jco) | 2 and is obtained from equation (1) as 

It 

| Vp(jw) j 2 = |G(ja>) | 2 1 v f (jw) | 2 (2) 


In the time domain, the autocorrelation of the particle velocity is 


V T) = 


=i im _L r » 

T -”2TL P 


(t)Vp(t + T)dt 


With zero time delay 
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y°> = t-“ ~ £ ''p<« dt 

which is just the definition of the mean square. So, 

yo> = <^> (a 

But, by the Wiener -Khintchine theorem, the autocorrelation is related to the energy 
spectrum by the inverse Fourier transform 



So, from equation (3), 



Similarly, the mean square of the fluid velocity is 

/ ’ O0 

|v f (jw) 1 2 du> (5) 

oo 

Substituting equation (2) into (4) : 



Hence, 
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J | G(ju>) | 2 1 v f (ja>) j 2 dc*j 

<v f 2 > 

f (v f (ju>) 1 2 dco 

y-O o 



Equation (6) is a weighting of the transfer function to the spectral distribution of 
energy in the flow; that is, the numerator of equation (6) is a measure of the particle’s 
ability to track the flow (through the transfer function) with consideration given to the 
actual distribution of frequencies present in the flow (through the energy spectrum of the 
flow). For particle trackability through a frequency bandwidth a> u , 



The energy spectrum of the main flow can be either measured directly or calculated 
from a measured autocorrelation by way of the Fourier transform 



Hence, to obtain the relative trackability of the particle, it is required to find the cor- 
rect transfer function G(jw). To determine the transfer function, the governing equation 
for the particle in the flow must be solved in the frequency domain. 
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Particle Motion in General Fluid Flow 


The general equation for the motion of a small spherical particle in a fluid, derived 
by Basset, Oseen, and Boussinesq (ref. 2), is 



where 

mass of spherical particle of diameter d with density p 

mass of equivalent spherical fluid particle with density p f 
dynamic viscosity of fluid 
viscous drag force 

external potential force (the weight due to gravity minus the upward buoyant 
force) 

The four basic assumptions in equation (8) are 

(1) The particle is spherical. 

(2) The particle is small when compared to the smallest wavelength of fluid motion. 

(3) The flow field is not perturbed by the presence of the solid particle. 

(4) Throughout the motion of the solid particle, the same fluid element remains in 
its neighborhood; that is, the pathlines of the solid particle and the fluid coincide. 

For a 1-micrometer particle stationary in an airstream at 1 atmosphere (14.7 psia) 
and 27° C (80° F) with a mean velocity of 30 meters per second (98.4 ft/sec), the maxi- 
mum particle Reynolds number is 



Re = III _ (30 m/sec) (lxlQ~^ m) 
1 ’ 1 . 57xl O' 5 m 2 /sec 


= 1.91 


Strictly speaking, the Stokes flow regime, where drag force is proportional to velocity 
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rather than to velocity squared, is valid only to a Reynolds number of about 1. How- 
ever, even at a Reynolds number of 5, the deviation in drag coefficient from Stokes law 
is not great; and it can be assumed that the drag on the particle is proportional to the 
relative velocity between the particle and the fluid. The drag force, then, is given by 
Stokes law as 


D = 3jrpd(v f - v ) 

The last term on the right side of equation (8) is any external potential force. In 
most applications this will simply be the weight due to gravity minus the upward buoyant 
force. In this case, the motion of the particle in steady state due to this force will be 
constant and equal to its terminal settling velocity in the given fluid. Because of the 
linearity of equation (8), this vertical component of the particle motion will not be af- 
fected by the general fluid flow field, and the total motion of the particle will be a linear 
superposition of its vertical settling velocity and its general motion induced by the flow 
field. The last term on the right side of equation (8), therefore, can be dropped without 
any loss of generality . The resulting equation is 


~ Pp 4 p = 3 "Pf d(v f - V + v Of*! * l X PA - V 



The meaning of each term in equation (9) is as follows: The single term on the left 
side is the mass of the particle times its acceleration; that is, it is the force required 
to accelerate the particle. 

The first term on the right side of equation (9) is the viscous drag on the particle, 
assuming Stokes flow. The second term is the force on the particle due to the pressure 
gradient in the fluid. This pressure gradient is caused by the acceleration of the fluid. 
And the resulting force is, in a true sense, a buoyant force in a direction opposing ac- 
celeration; that is, in so far as the particle is concerned, the effect of an accelerating 
fluid field is the same as a gravity field. The third term is the force required to accel- 
erate the apparent mass of the particle relative to the fluid. The last term, derived by 
Basset, accounts for the deviation of the flow from steady state. 

Rewriting equation (9) gives 
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(10) 


where 


and 



Simplification of Flow Equation 

Equation (9) suggests that when Pp/pf the second, third, and fourth terms of 
the right side of equation (9) become negligibly small. The equivalent limit in equa- 
tion (10) is 0—0 with o?0 simultaneously remaining finite. This would leave only the 
first term on the right side of equation (10); that is, the only driving force would be a 
linear drag force. It is worth investigating the specific characteristics of the flow that 
would be required for this assumption to be valid. To do this, equation (10) must be 
nondimensionalized. As a reference velocity, let 

v ref = “> 

where u> is some frequency associated with the main flow. As a reference time, let 

‘ref ’ (5i) ' 1 

The nondime nsional velocities and time then become 
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and 


So, 


Also, 




v w * 

= uxiVj 



Therefore , 



u> 


Substituting equations (11) 


and (12) into equation (10) gives 


3 a d ^ + a0SW„ = apwdVj + 
de 


-2 



So^V 72 


Rearranging gives 


8 



(13) 


d0 


(9 ’ T9 


= -2-V f + /3 



It is now clear when the last two terms on the right side of equation (10) can be 
neglected and only the linear drag force retained. In the limit as 0 — 0 (i.e. , P^/Pf — °°) 
and at the same time the dimensionless frequency (co/at) becomes very small, the second 
and third terms on the right side of equation (13) become negligible in comparison with 
the remaining terms. More specifically, if (aj/of/3) = ^(1) or less as 0 — 0, the pre- 
dominant driving force is the linear drag force. As the flow departs from this condi- 
tion, the next to become signficant is the Basset force, and finally the drag force pro- 
portional to acceleration. In other words, for low-frequency flow with a large particle - 
to fluid-density ratio, only the linear viscous drag force need be considered. 

Returning to equation (10) with only the viscous drag force retained, 


v p + <*/3v p = a/3v f 


(14) 


Taking the Fourier transform of equation (14) yields 


jwvp(jw) + a$Vp(ia>) = a$Vf(jo>) 


So, 


~ a09 f ( ju) 

v a/3 + ja> 


The complex transfer function is 


G-Qw) = -2 = 

v f (ja>) oi(3 + jw 


aj3 _ (oigr - 

(oi0) 2 + w 2 
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Thus, 



and from equation (6) , 



For the assumed model, as given by equation (14), equation (16) is valid regardless 
of the nature of the fluid flow field: turbulent or laminar, time varying or constant. 

The simplest case is that of a sinusoidally varying flow field 

v f (t) = v f cos w Q t (17) 

The energy spectrum of the flow field is given by the Fourier transform of equation (17) 

«?"[v f (t)] =^[ v f cos w Q t] = 7 t[6(w - w Q ) + 6(w + u> 0 )] = v f (ja>) 

where 6 is the Dirac -delta function. So, from equation (16), 
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+ 0 ) 


2 

o 


GjOco) 


(18) 


Therefore, for a simple sinusoidal velocity field, the particle trackability , or the true 
frequency response, is the magnitude of the transfer function. Although no real flow 
will be described by an equation like (17), the magnitude of the transfer function as a 
function of frequency does serve somewhat as a measure of relative particle trackability 
under different flow conditions, in the absence of any spectral knowledge of the fluid 
flow. The quantity 1 /(o?/3) can be considered as a crude time constant: the smaller . 
1/{oi0), the better the particles will track for given fluid flow conditions. 

A plot of equation (18) as a function of (u) o /af3) is shown in figure 1. To more 
graphically illustrate the effect on trackability of particle size and frequency, equa- 
tion (18) is plotted in figure 2 for the specific case of aluminum oxide (A^Og) particles 

o o e 9 

(pp = 3.5 g/cm ) in an airstream (p^ = 1.2 kg/m , = 1.85x10 N-sec/ncO with sinu- 

soidally varying velocity, for various particle diameters. From the figure, it can be 
seen that there is a severe drop in particle trackability with frequency. For example, 
the velocity of a 0. 5-micrometer particle will track within 1.5 percent of the true fluc- 
tuating velocity of the fluid at 10 kilohertz, with correspondingly greater errors for 
larger particles and higher frequencies. 

More realistically, the fluid flow is not oscillating at a single frequency but pre- 
sents a broad spectrum of frequencies for the particle to follow. It is not important that 
the particle can track the fluctuating velocity of the fluid within 20 percent at 50 kilo- 
hertz if, for example, most of the turbulent energy occurs at frequencies below 10 kilo- 
hertz. 

The effect of this weighting of the transfer function with the energy spectrum of the 
flow is best shown by example. Laurence (ref. 3) has shown that the energy spectrum 
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of free jet turbulence can be given by an expression of the form 


| v-(jco) 


1 2 _ (1 + B)tt 


(If 

(ifl Nij 


where w is some characteristic frequency and B is a constant. Both to and B can 
be obtained from the autocorrelation of the flow at the point of interest. Specifically, 
oj and B are obtained from the equations 


. -Ba+B>/U-B> 


i. lnB 

r 0 (B - 1) 


(20b) 


where R f (r) min is the minimum value of the autocorrelation and r Q is the value of the 
time delay at the first zero crossing of the autocorrelation. 

Putting equation (19) into (7) and using equation (15) yield 
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Equation (21) was evaluated numerically for the specific case of particles 

(p = 3. 5 g/cm ) in an airstream (p f = 1.2 kg/m , = 1 . 85x10 N-sec/m ) using the 

Laurence data for a point 1.14 diameters downstream of the nozzle exit plane and 1/2 
nozzle diameter from the jet axis. 

The result is plotted in figure 3. Some caution should be exercised in the inter- 
pretation of figure 3. The figure does not indicate, for example, that a 1 -micrometer 
particle will faithfully reproduce the rms velocity of this flow within 2 ^ percent at 
1 megahertz. It does indicate, however, that for this particular flow, if information is 



Figure 3. - Accuracy of reproduction of rms velocity in a 8.9-centimeter (3. 5-in. ) diameter free jet for AI 2 O 3 par- 
ticles of various diameters. (Data from ref. 3. ) 


acquired through the entire spectrum to 1 megahertz, the measured rms velocity of the 
particle will be within 2j percent of the true rms of the fluid in that bandwidth. The 
reason for this is simply that the energy spectrum for this point, as given by equa- 
tion (19) , contains little significant energy beyond the breakpoint of the corresponding 
transfer function given by equation (15); that is, the energy spectrum drops more rap- 
idly with frequency than the transfer function. 
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Comparison with Complete Solution 


Within the limits of the assumptions, equation (21) is a fair representation of the 
trackability of a single particle in a turbulent flow. For best trackability , however, it 
is obvious that the particle density should be as close as possible to the fluid density. 
For example, the same AlgOg particles in a water flow give a Pp/p£ of 3.5. Under 
such circumstances, the assumption of p p /p f » 1 would no longer be valid, and the 
full equation (10) would have to be solved to determine trackability. 

The solution to equation (10) for a stationary flow of infinite extent (i.e. , no solid 
boundaries) in terms of the transfer function, as presented in reference 2, is 


O' + 


G 2 (ju>) = -t 


a + 



and 



and 



As the flow becomes very slowly varying and approaches dc (i.e. , u>/a —0), then 
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|G 2 (jo>)j -1. Further, as 0 - 1 (i.e., p p /p f - 1), then |G 2 (jco)j -1 for all values of 
u/ci. In other words, for slowly varying flows or for neutral density particles, perfect 
trackability is approached. However, if p p /p £ <l,then 0 > 1 and |G 2 (jco)[ > 1. So, 
under these circumstances the particle motion overshoots the fluid motion. Hence, at- 
tempting LDV measurements with very light particles (relative to the fluid) can intro- 
duce errors just as large as with relatively heavy particles. 

That the solution to equation (10) approaches the solution to equation (14) when the 
ratio of particle density to fluid density gets very large can be seen by applying this 
limit to the respective transfer functions. First, rewrite equation (22) as 



M |G 2 (1 “>I = 

finite 

a0 



which is the same as equation (15). 


The solution, then, can be written as 
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( 23 ) 


Equation (23) is an exact solution to equation (10); and either equation (23) can be 
numerically evaluated if an expression for v^(jw) is available, or the transfer function 

alone can be used as a rough measure of trackability . When ^(jcu) is 
available, though , the numerical computation of equation (23) is somewhat lengthy even 
with a high-speed computer. For an accurate error analysis, equation (23) would have 
to be evaluated at each point of interest in the flow under consideration, and many such 
computations require a large amount of computing time. The same computation using 
the transfer function for the simplified analysis (i.e. , eq. (16)) is considerably faster 
and it is therefore of some interest to determine when equation (16) can be used with 
accuracy and when equation (23) must be used. 

The analysis of the nondimensional equation (13) indicated that for pp/p f -°°, the 
reduced equation (14) could be used if 0(u/ct) = ^(/3) (i.e. , at low frequency). A plot 



Figure 4. - Variation of transfer function ratio with dimensionless frequency for various particle- to fluid-density 
ratios. 
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of the ratio of the simplified transfer function (eq. (15)) to the full transfer function 
(eq. (22)) against <x>/ar/3 is shown in figure 4 for various values of pp/pf The value 
of oo/afi at which a given curve deviates significantly from unity represents the fre- 
quency above which the full equation (23) must be used. These curves, however, must 
always be used in conjunction with some spectral knowledge of the fluid flow. For a 
given fluid and seeding material, a curve in figure 4 might indicate that the simplified 
transfer function can be used to 100 kilohertz. Nevertheless, if the general fluid flow 
contains significant energy above 100 kilohertz, equation (23) must be used. 


CONCLUDING REMARKS 

It was shown that for the limiting case of large particle- to fluid-density ratio and 
low dimensionless frequency, the reduced flow equation can be used to compute particle 
trackability . The question then arises as to why one would use particles whose density 
is large compared to the fluid it must track, since the larger the deviation from neutral 
density, the larger the tracking error. One answer is that there is not always complete 
freedom in the choice of seeding particles. For example, many (if not most) applica- 
tions would require seeding materials to be nontoxic, noncorrosive, nonabrasive, etc. 
This would eliminate a large number of smokes and similar aerosols. Similarly, if the 
flow environment is severe, such as high temperature, then many of the synthetically 
produced materials, such as polystyrene particles are eliminated. Finally, when cost, 
availability, and ease of flow seeding are considered, the number of practical seeding 
materials is greatly reduced and, for a given application, only a relatively dense seed- 
ing material might be suitable. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 26, 1972, 

132-80. 



APPENDIX - SYMBOLS 


B 

D 

d 

f 

G(jaj) 

Gj(jw) 

G 2 (jw) 

j 

C 

P g 

H(r) 

Re 

T 

t 

t’ 

V 
v 
v 

V 

a 

P 

6 

6 » 

6’ 

X 

P 

v 


constant defined by eq. (20a) 

drag force 

particle diameter 

Fourier transform 

frequency, Hz 

complex transfer function 

complex transfer function for eq. (15) 

complex transfer function for eq. (11) 

fl 

order of magnitude 
force due to gravity 
autocorrelation function 
Reynolds number 
limit of integration, sec 
time, sec 

dummy variable of integration, eq. (8) 

dimensionless velocity 

velocity 


dv/dt 

velocity in frequency domain 

2 

time constant, 12^/d 
density ratio, 3p f /|2p p + p^J 


Dirac -delta function 
dimensionless time 

dummy variable of integration, eq. (13) 
wavelength 
dynamic viscosity 
kinematic viscosity 
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p density 

r delay time, sec 

Tq delay time at zero crossing of autocorrelation function 

functions defined by eq. (22a) and (22b) 
oj circular frequency, rad/sec 

dS characteristic circular frequency defined by eq. (20b), rad/sec 

aT reference circular frequency, rad/sec 

Subscripts: 

f fluid 

p particle 

o initial time 

ref reference 

rms root-mean-square 

u upper limit 

Superscripts: 

* complex conjugate 

( ) time average 
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